January 12, 2013 12:36 



WSPC - Proceedings Trim Size: 9in x 6in ws-proc59x6 



1 



Nuclear matter, nuclei, and neutron stars in hadron and 
quark-hadron models 

S. Schramm 

Frankfurt Institute for Advanced Studies, Goethe University, 
Frankfurt, D-60438, Germany 

* E-mail: schramm&fias .uni-frankfurt. de 

V. Dexheimer 

Gettysburg College, 
300 North Washington Street Gettysburg, PA 17325 USA 
E-mail: vantoche@gettysburg.edu 

R. Negreiros 

Frankfurt Institute for Advanced Studies, Goethe University, 
Frankfurt, D-60438, Germany 

* E-mail: negreiros&fias.uni-frankfurt. de 

T. Schiirhoff 

Frankfurt Institute for Advanced Studies, Goethe University, 
Frankfurt, D-60438, Germany 
* E-mail: schuerhofJQth.physik. uni-frankfurt. de 



We develop a unified model of hadrons and quarks. Within this approach we 
investigate the phase structure of the model as function of temperature and 
chemical potential. Computing the equation of state of cold matter we deter- 
mine neutron and hybrid star masses and radii. In an extension of the inves- 
tigation we consider the cooling behavior of the compact stars and derive a 
general relation between the star's mass and rotation and its cooling behavior. 
Finally we study the effect of A resonances for star matter, especially with 
respect to possible solutions of stars with small radii. 
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1. Introduction 

A central topic in modern nuclear physics is the stiidy of the general behav- 
ior of strongly interacting matter under extreme conditions. One expects 
a transition from the effective degrees of freedom in vacuum and ground 
state matter, hadrons, to a phase of quarks and gluons at high tempera- 
tures and/or densities. This can be described by the restoration of chiral 
symmetry and the deconfinement transition. The actual structure of the 
transition is still unknown, perhaps with the exception in the case of van- 
ishing chemical potential, where lattice results show a cross-over transition. 
Some lattice results also point to a critical end-point of a first-order tran- 
sition line in temperature and chemical potential.^ In order to describe the 
matter, a main challenge for thciory is to develop models that can cover 
all the ranges of temperatures and densities. This requires a model that 
is able to reproduce both nuclear matter ground state properties and the 
quark and gluon phase beyond the phase transition. In the following we 
develop an effective theoretical description along these lines. Results for 
neutron star properties and their cooling behavior within this approach 
will be presented. In addition we study the influence of the A baryons on 
stellar properties. 

2. Effective Chiral Model 

In order to describe hadronic physics as well as quark degrees of freedom we 
couple an effective chiral flavor-SU(3) model for hadrons^ with a modified 
version of the so-called PNJL quark model. The model includes the lowest 
SU(3) multiplets of hadrons and quarks. Baryonic and constitutent quark 
vacuum masses are generated via spontaneous symmetry breaking by the 
coupling of the particles to the scalar fields, which attain non-zero vacuum 
expectation values. Explicitly the Lagrangian reads:^ 

L = Lxin + Lint + Lself + LsB - U, (1) 

Here, the interaction term between baryons, quarks and mesons is given by 

Lint = - X! M'yoidiujt^ + 9i<j><l> + gipT3p) + M*]'ipi, 

i 

with the effective baryon and quark mass 

M* = 9ia<y + 9iST3S + Qi^C + SMi . (2) 

SMi is a small explicit mass term. The self interactions of the scalar fields, 
generating the spontaneous symmetry breaking, and of the vector fields are 
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contained in 

Lseif = -limluj'^ + mlp^ + mlcf)'^) 
+fco(tT2 + + 52) + fci(a2 + + 

+k2 + f + 3C72<52 + + ^3(^2 _ ^2)^ 

+fc4 ln^4#^. (3) 
Finally, the explicit symmetry breaking term Lsb lias the form 

LsB = mlUa + (V2mlfk - /-) C • (4) 

The mesons included are the vector isoscalar fields oj and (f), the vector 
isovcictor p, the scalar isoscalars a and (non-strange and strange quark- 
antiquark state, respectively) and the scalar-isovcctor state 6. The isovector 
mesons affect isospin asymmetric matter, which makes them relevant for the 
study of compact stars. The various couplings are fitted to reproduce the 
vacuum masses of the baryons and mesons, nuclear saturation properties 
(density po = 0.15 fm~'^, binding energy per nucleon B/A = —16.00 MeV, 
nucleon effective mass = 0, 67 Mn, compressibility K = 297.32 MeV), 
asymmetry energy {Esym = 32.50 MeV), and reasonable values for the 
hyperon potentials {Ua = -28.00 McV, = 5.35 MeV, Us = -18.36 
MeV). The numerical values of the couplings can be found in Ref. 2-4. 

In order to be able to study the transition from hadronic to quark de- 
grees of freedom the model incorporates a Polyakov loop field $ in analogy 
to the PNL models for quarks. ^'^ The field has an effective potential of the 
form 

+a3r4log(l-6$2 + 8$3_3$4)^ 

depending on the temperature T and the baryochemical potential p. The p 
independent terms have been fitted, so that thermodynamical quantities of 
quenched lattice calculations agree well with the model results over a wide 
range of temperature. The other two parameters of the potential have been 
adjusted in order to obtain a phase diagram that has a first-order phase 
transitions at high densities which ends in a critical end-point at a tem- 
perature T and chemical potential p in agreement with results from lattice 
calculations.^ In order to suppress quarks at low densities and temperatures 
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and hadrons in the deconfined phase an efffective mass shift is introduced 
of the form:^ 

5M*s = 9B^^^ , SM* = gg^{l - $) . (6) 

3. Results 

In an earher investigation, where we fine-tuned the couplings to properties 
of finite nuclei, we computed nuclei in a two-dimensional calculation includ- 
ing nuclear deformation^ . The corresponding equations of motion for the 
mesonic fields as well as the Dirac equation for the nucleons were solved on 
a spatial grid. Computing all nuclei with known binding energies an overall 
accuracy of the energies could be achieved with an average error for nuclei 
with 50 or more nucleons of 

eB(A > 50) w 0.21% . (7) 

This result compares very wcill with standard rclativistic and non- 
relativistic mean-field nuclear structure models. As can be seen in Fig. 1 
magic numbers can be reproduced as well as the measured deformations of 
many nuclei. Thus this shows that a. more ffeneral hadronic model approach 
can be used t 




Fig. 1. Nuclear chart as calculated in the effective cliiral model. The gray level encodes 
the deformation (absolute value) of the corresponding nucleus. One can clearly observe 
lines of spherical nuclei corresponding to magic numbers. 



Investigating the full model and tuning the couplings in order to repro- 
duce high-mass neutron stars as well as having reasonable nuclear ground 
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state values^ we minimize the grand canonical potential in mean-field ap- 
proximation following from the Lagrangian discussed in the previous sec- 
tion.'^ Fi gure 2 shows the resulting phase transition curve for symmetric 
nuclear matter and star matter in j3 equilibrium. The critical end-point oc- 
curs at Tc ~ 167 MeV and /ic ~ 354 MeV in agreement with Ref. 1. With 
our choice of parameters the phase transition to quarks takes place at 4 
times ground state density for vanishing temperature. Since the model fea- 
tures a realistic nuclear ground state the figure also shows a quantitatively 
reasonable first-order liquid gas phase transition. 
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Fig. 2. Phase diagram as function of temperature and baryon chemical potential. The 
first order phase transition stops at a critical end point as suggested in Ref. 1. At T = 
the transition takes place at 4 times nuclear matter saturation density. The phase 
transition for symmetric and star matter is shown. 

Computing the resulting field values one can determine energy density 
and pressure for stellar matter. Using the equation of state we solve the 
Tolman-Oppenheimer-Volkov equations for non-rotating, spherical neutron 
stars. The resulting masses and radii of the stars are shown in Fig. 3. 
Excluding quarks the maximum attainable mass is about 2.1 solar masses. 
Taking quarks into account the maximum mass drops to « 2 solar masses. 
Note that, when quarks are populated, the region of stable stars is cut 
off as the equation of state in this regime becomes softer and cannot sup- 
port stable high-mass star solutions any more. This is true for a Maxwell 
construction of the phase transtion assuming local charge neutrality of the 
matter. In a Gibbs construction, which only assumes global charge neu- 
trality, allowing for a mixture of quark and baryonic phases with opposite 
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charge, this sharp cut off is softened as can be seen in the inset of the figure. 
As a result at the highest star masses a core of about 2 km size develops 
that contains a mixture of quarks and baryons. This can also be observed in 
Fig. 4, which shows the baryon number densities of baryons and quarks as 
function of chemical potential. Around the phase transition, the plot shows 
this mixture of etafps 
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Fig. 3. Masses and radii of neutron stars. The results for purely baryonic stars and for 
stars including the quark phase are shown. The inset illustrates the difference between 
a Maxwell and a Gibbs construction of the hadron-quark phase transition. 
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Fig. 4. Baryon number densities of particles for star matter as function of baryon chem- 
ical potential. The mixture of baryons and quarks results from the Gibbs construction of 
the first order phase transition assuming a mixed phase with global charge conservation. 
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One of the most interesting observables concerning neutron stars is their 
temperature compared to an age estimate of the star, defining the coohng 
behavior of the neutron star. To study this feature we make use of the 
equation of state of the hadron-quark model for stellar matter as described 
before. Incorporating this equation of state in a cooling simulation of the 
star we compute the cooling behavior of stars with different masses. Fig- 
ure 5 shows the result of the calculation. Compared to the measurements 
also shown in the figure, one can conclude that a good agreement with data 
can be achieved if the corresponding (unknown) star masses are roughly 
below 1.2 solar masses. For higher masses the model predicts much faster 
cooling rates. The main reason for this behavior is the on-set of direct Urea 
processes for heavy stars. The most important processes are 

n ^ p + e'^ + , p + ^ n + Ve ■ (8) 

These processes generate fast cooling via neutrino emission. In a more gen- 
eral analysis we take into account rotation of the star. As rotation reduces 
the central energy density of the star also the threshold of the direct Urea 
process shifts to higher masses. The result of this analysis is shown in Fig. 
6. Here one can read off the ranges of frequencies and masses of the stars 
that allow for the direct Urea process to take place, which leads to fast 
cooling of the star. 




Fig. 5. Cooling behavior as function of time for compact stars as calculated in this 
model. Observational data are also shown. One can observe a good agreement with the 
data assuming low star masses below 1.2 solar masses (assuming non-rotating stars). For 
higher masses the direct URCA process leads to enhanced cooling. 
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Fig. 6. The figure illustrates the range of star masses and rotation frequencies for which 
the Urea process takes place (DU process allowed). The shaded region at high frequencies 
is not physically realized as the stars would be instable with respect to mass shedding. 



In a separate investigation concentrating on a purely baryonic star, for 
simplicity, we investigate the possibility of the existence of neutron stars 
with rather small radii around 10 km and below, as has been suggested in 
an analysis of masses and radii of three neutron stars discussed in Ref. 11. 
Here we also include the spin-3/2 baryonic multiplet as degrees of freedom 
as has already been studied in Ref. 3. As with other baryons the couplings 
of the particle to the scalar fields are fixed by reproducing baryonic vaccum 
masses. The vector interaction strength is less clear. Using vector domi- 
nance arguments a natural choice is to take the same value as in the case of 
the corresponding baryonic octet state with the same number of light and 
strange quarks. Starting from this assumption we allow for a small reduc- 
tion from this canonical value. The resulting star masses and radii are 
shown in Figure 7. One can see that for a moderate change in the vector 
strength of about 10 percent {ry = 0.9) one obtains good agreement with 
the values obtained in the analysis of the observational data. Thus it is 
possible to reproduce stars with relatively small radii in a purely baryonic 
model without introducing a more exotic solution like a quark star. The 
population densities of the various baryons as function of chemical poten- 
tial are shown in Fig. 8. It is interesting to observe a very different structure 
compared to calculations not including the A resonances. In fact, beyond 
4 times nuclear ground state densities the stellar matter is dominated by 
A's. 
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Fig. 7. Neutron star masses and corresponding radii. The strength of the vector- 
interaction of the Delta resonance is varied. Also shown are the results of an analysis 
of three observed neutron stars. denotes the value of the vector coupling of the A 
relative to the nucleonic value, i.e. r„ = gAu/SJVij. The model calculation can reproduce 
the observational data for a relative vector strength = 0.9. 




Fig. 8. Particle number densities as function of density including A resonances. One 
can observe that the As dominate other particles at densities p > 4po. 



4. Summary 

We developed a unified hadron-quark model with the correct degrees of 
freedom on both sides of the QCD phase transition. Within this approach 
we were able to reproduce a phase transition line with a critical end point for 
values of temperature and chemical potential that were predicted by some 
lattice gauge calculations. Using these model parameters we studied neutron 
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star properties and found heavy stars with maximum mass vahics of around 
2 solar masses. The core of the star contained a mixture of baryons and 
quarks. Including a study of the cooling of the star we observed agreement 
with data for low-mass stars. In the case of high masses the direct Urea 
process led to fast cooling rates. Wo determine the general parameter region 
of stars with respect to rotational frequency and mass that did not allow 
for the Urea process to take place suppressing fast cooling of the stars. 

In an investigation of purely baryonic stellar matter but taking into ac- 
count the baryonic spin-3/2 multiplet, in particular the A resonances, we 
showed that one is able to reproduce stars with small radii, which were sug- 
gested by the analysis of observational data. This was achieved by a slight 
tuning of the vector coupling strength of the spin-3/2 particles without the 
need for exotic particles. 
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